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ABSTRACT 

Lindal et al. (1987, J. Geophys. Res. 92, 14987-15001) presented a range of temperature and 
methane profiles for Uranus that were consistent with 1986 Voyager radio occultation measurements 
of refractivity versus altitude. A localized refractivity slope variation near 1.2 bars was interpreted to 
be the result of a condensed methane cloud layer. However, models fit to near-IR spectra found particle 
concentrations much deeper in the atmosphere, in the 1.5-3 bar range (Sromovsky et al. 2006, Icarus 
182, 577-593, Sromovsky and Fry 2008, Icarus 193, 211-229, Irwin et al. 2010, Icarus 208, 913-926), 
and a recent analysis of STIS spectra argued for a model in which aerosol particles formed diffusely 
distributed hazes, with no compact condensation layer (Karkoschka and Tomasko 2009, Icarus 202, 
287-309). To try to reconcile these results, we reanalyzed the occultation observations with the He 
volume mixing ratio reduced from 0.15 to 0.116, which is near the edge of the 0.033 uncertainty range 
given by Conrath et al. (1987, J. Geophys. Res., 15003-10). This allowed us to obtain saturated 
mixing ratios within the putative cloud layer and to reach above-cloud and deep methane mixing 
ratios compatible with STIS spectral constraints. Using a 5-layer vertical aerosol model with two 
compact cloud layers in the 1-3 bar region, we find that the best fit pressure for the upper layer is 
virtually identical to the pressure range inferred from the occultation analysis for a methane mixing 
ratio near 4% at 5° S. This strongly argues that Uranus does indeed have a compact methane cloud 
layer. In addition, our cloud model can fit the latitudinal variations in spectra between 30° S and 20° 

N, using the same profiles of temperature and methane mixing ratio. But closer to the pole, the model 
fails to provide accurate fits without introducing an increasingly strong upper tropospheric depletion 
of methane at increased latitudes, in rough agreement with the trend identified by Karkoschka and 
Tomasko (2009, Icarus 202, 287-309). 

Subject headings: Uranus, Uranus Atmosphere; Atmospheres, composition, Atmospheres, structure 


1. INTRODUCTION 

The existence of a thin methane ice cloud in t he at - 
mosphere of Uranus was inferred bv lLindal et al.1 ( 1987 ), 
henceforth referred to as L87, from their analysis of Voy¬ 
ager 2 radio occultation measurements. They also de¬ 
rived a suite of temperature and methane profiles that 
were all consistent with their measurements, including 
Model F, which had the greatest deep methane mixing 
ratio (4%), and their preferred Model D, which had a 
deep mixing ratio of 2.3% and an above-cloud relative 
humidity of 30%. (Here we use the term relative humid¬ 
ity of methane to refer to the ratio of its partial pressure 
to its saturation pressure at the same temperature, and 
the mixing ratio referred to is the volume mixing ratio or 
VMR.) A methane cloud layer near the 1.2-bar level in¬ 
ferred by L87 has been used successfully in the analysis of 
o bservatio ns in the visible spectral range. For example, 
iRaees et al.1 (|1991l ) incorporated a methane cloud layer 
into their models of Voyager 2 imaging observations, find¬ 
ing modest optical depths of 0.66±0.18 at 22° S (assumed 
to be independent of wavelength), and about three times 
that level at 65 ° S, while assuming a deep methane 

1 University of Wisconsin - Madison, Madison WI 53706 
t Based in part on observations with the NASA/ESA Hub¬ 
ble Space Telescope obtained at the Space Telescope Science 
Institute, which is operated by the Association of Universities 
for Research in Astronomy, Incorporated under NASA Contract 
NAS5-26555. 


mixing ratio of 4% (consistent with L87 Model F). On 
the other hand, an analysis of hydroge n S(0) and S(l) 
quadrupole line features bv I Baines et al.1 (Il995h . which 
also incorporated a methane cloud near 1.2 bars, found 
its opacity (weighted to high latitudes) to be about 0.4 
at 0.6 /im, while inferring a deep methane mixing ratio 
of 1.6%. Neither of these authors tried to constrain the 
pressure of the cloud, however, so that consistency with 
occultation results was not fully established. 

Recently, a serious challenge to the exis¬ 
tence of the meth ane cloud layer was made by 
1 Karkoschka and Tomaskol ( 2009 )1. henceforth referred to 
as KT2009, based on their analysis of spatially resolved 
0.3-1 fim spectra obtained from 2002 observations by 
the Hubble Space Telescope Imaging Spectrograph 
(STIS). They concluded that the most significant cloud 
opacity concentration was in a layer from 1.2-2 bars, 
with particles uniformly mixed with the gas in this 
layer, which had wavelength-independent optical depths 
between 1.2 and 2.2. They argued for no localized CH 4 
condensation layer at all, but instead for the existence 
of a global thick and diffuse tropospheric haze similar 
to that observed on Titan. This seemed to confirm the 
analysis of near-IR spectral observations, which had 
already questioned the existence of a methane ice cloud 
near 1.2 bars. 

From an analysis of the iFink and Larsonl (11979D near- 
IR spectrum, which made use of improved methane ab- 
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sorpti on coefficients (llrwin et a l. 2006). [Sromovsky et al.1 
(2006) obtained a cloud layer near 2 bars for the L87 
Model D profiles and at 1.5-1.7 bars for their Model F 
profiles. A subseq uent a nal ysis o f 2004 near-IR imag¬ 
ing observations bv lSromovskv and Frvl ( 2007 ) concluded 
that the methane relative humidity should be near 60%- 
100% above the nominal cloud region (Models D and 
F had 30% and 53% respectively), and at low latitudes 
found no need for a cloud layer in the 1.2-2 bar region. 
Further refineme nts o f metha ne absorpt ion models for 
the near-IR (jKarkoschka and Tomaskoll2010D did not en¬ 
tirely fix these discrepancies! llrwin et al.l (j2010( l used 
these improved methane coefficients and the L87 Model 
D T(P) profile and above-cloud methane mixing ratio, 
but used 1.6% instead of 2.26% below the putative cloud 
layer. With these assumptions they obtained a low- 
latitude cloud density peak near 2.5 bars. Even using 
the L87 Mo del F T(P) profile and a 4% deep methane 
VMR, llrwin et all (1201 Of ) still obtained a cloud peak that 
was too deep to reach the methane condensation level, 
though the cloud pressure was then elevated to about the 
1.7-bar level, which puts the peak in the middle of the 
main aerosol layer of KT2009. 

More methane seems to be required to bring the cloud 
pressure inferred from spectral observations to the same 
level as inferred from the refractivity profile. However, 
according to L87, their Model F is an upper limit on 
methane amounts, and they argued that Model D is re¬ 
ally preferable. That solution has a deep methane mix¬ 
ing ratio of 2.3% by volume, a cloud layer between 1.15 
and 1.27 bars, and an above-cloud methane humidity of 
30%. Other profiles that satisfy the occultation mea¬ 
surements have deep methane mixing ratios varying from 
zero to 4%, and above-cloud humidities varying from zero 
to 53%, while the in-cloud humidities vary from zero 
to 78%. Model D was preferred for three reasons: (1) 
it provides the best agreement with IRIS observations 
sampling the above-cloud region, (2) it has the high¬ 
est in-cloud humidity levels, (3) it yields a deep mix¬ 
ing ratio (2 .3 %) th at is in close agreement with that of 
lOrton et al.l ( 1986 ). The first rea son is weak because 
the IRIS observations in question dConrath et al.1 fl987h 
are at a large zenith angle and rather uncertain. The 
second is weakened by the fact that only the nominal he¬ 
lium volume mixing ratio was considered, and that can 
strongly affect the methane humidity, as we will show 
here. The third reason is questionable because the Orton 
et al. temperature profile derived assuming 2% methane 
does not agree with the occultation profile using vir¬ 
tually the same mixing ratio. Much stronger external 
constraints are available from spectral observations in 
the CCD (~0.3 - 1 (im) wavelength region, as shown by 
KT2009 and by the analysis presented here. In prior use 
of these spectral constraints however, there has been an 
unjustified deviation from occultation constraints in both 
thermal and methane profiles, in which a thermal profile 
derived for one me thane profile is used for a very dif¬ 
ferent mixing ratio (llrwin et al.ll2010t iBaines et al.l ll995: 


iSromovskv and Frvl 2007t) amlabove-cloud methane pro¬ 

files have been used that e xceed a ll of the occu ltation so¬ 
lutions, e.g. KT2009 and IBaines et al.1 (1199511 . KT2009 
also inferred that the methane mixing ratio varies with 
latitude, which raises additional questions about the ef¬ 
fects of corresponding density variations with latitude. 


To summarize, where spectral observations have been 
used to test the location of cloud layers on Uranus, the 
inferred locations are considerably deeper than implied 
by the occultation observations. And the spectral con¬ 
straints on the methane mixing ratio range from a low 
of 1.6±^% bv IBaines et al.1 (119951) to 4% bv lRages et al.1 
dHH, and include latitude dependent values between 
these values inferred by KT2009. Further, with the ex¬ 
ception of Rages et al., prior modelers have not followed 
the occultation constraints on the the vertical distribu¬ 
tion of methane. This motivates our efforts to redo the 
occultation analysis with consideration of a wider range 
of solutions, and to examine more carefully the plausi¬ 
bility of a compact methane cloud layer on Uranus. 

In the following, we pursue the point of view that 
the occultation measurements of sudden slope changes 
in refractivity do indicate a region of sudden changes in 
methane mixing ratio, which are indicative of the con¬ 
densation level of methane, and possibly of a thin region 
containing cloud particles. We describe a reanalysis of 
the occultation measurements that can actually achieve 
saturated vapor pressures in the same region as the sud¬ 
den changes in refractivity slope. We also find solutions 
with high methane amounts at and above the cloud level 
that are consistent with the adopted profile of KT2009. 
After finding a range of solutions with the desired char¬ 
acteristics, we then constrain these solutions by calcu¬ 
lating spectra for compact cloud layer models and com¬ 
paring the pressures inferred from matching spectra to 
the pressures inferred from the occultation analysis. We 
conclude that a methane cloud layer at the occultation 
pressure is consistent with the spectral observations, but 
that most of the cloud opacity is concentrated in a deeper 
layer that was not detected by the occultation measure¬ 
ments. We also confirm the conclusions of KT2009 that 
the methane is strongly depleted at high latitudes, but 
to shallow depths. 

2. APPROACH TO OCCULTATION ANALYSIS 

We first describe the basic methods of analysis, how we 
reconstruct the refractivity profile, then use that profile 
to validate our methane retrieval methods by comparison 
with results of L87. 

2 . 1 . Physical Basis and Methods of Occultation Analysis 

The occultation measurements, after accounting for 
observing geometry, can be reduced to refractivity as a 
function of altitude, where refractivity is defined to be 
N = (index of refraction -l)xl0 6 . If the molecular com¬ 
position is known, then refractivity can be converted to 
number density and mass density. The pressure can then 
be determined by integrating the product of mass den¬ 
sity and gravity, assuming hydrostatic equilibrium. From 
pressure and density, the equation of state of the gas can 
then be used to infer temperature. The main complica¬ 
tion is that the variable distribution of methane is not 
known and cannot be directly inferred from the observa¬ 
tions. This allows a range of T(P) solutions that depend 
on what is assumed about the methane distribution. 

The procedure followed by L87 was to first select a 
molecular weight that yielded a thermal profile near the 
tropopause that was consistent with IRIS thermal in¬ 
frared observations. The only significant constituents at 
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that level of the atmosphere are hydrogen and helium, so 
that the molecular weight was determ ined by their ratio, 
which was taken to be the iConrath et all (Il987ll value of 
He/H 2 =15/85. L87 did not consider other He/H 2 ratios, 
even though the quoted uncertainty is large enough to 
permit substantially different profiles, as will be shown 
in Sec. [3] The next step in the procedure was to select 
altitudes bounding the cloud layer. These altitudes were 
not stated by L87, but are approximately between 5 and 
7 km below the 1 bar level determined for the D model. 
These are the rough locations of the rapid changes in re- 
fractivity slope. While the altitude relative to the center 
of the planet is fixed for all the profile solutions, the pres¬ 
sure varies somewhat from one solution to the next, so 
that the altitude above the 1 bar level also varies slightly. 

To constrain the methane profile L87 generally as¬ 
sumed a constant relative methane humidity above the 
cloud layer and used the tropopause mixing ratio to 
set the constant stratospheric mixing ratio. Within the 
cloud region, L87 state that the temperature lapse rate 
was set equal to the wet adiabatic lapse rate, and ad¬ 
justed the number density and temperature to match the 
refractivity profile. However, the model D T(P) profile of 
L87 does not match the wet adiabat within the assumed 
cloud layer. Instead, the profile matches a weighted av¬ 
erage of the form 

(< dT/dz) cM = (< 4T/dz) dry x (1 RH) + (1) 

(dT/dz) wet X RH 

where RH is the relative humidity. This weighted av¬ 
erage is consistent with the suggestion that the occul- 
tation sampled a broad horizontal region in which the 
average humidity was less than expected for a uniform 
cloud layer. That is also offered as an explanation for the 
sub-saturated humidity levels that were obtained in the 
cloud layer. The relatively smooth I/F profiles observed 
on U ranus (|Sromovskv and Frvl 120071 : iSromovskv et al.1 
l2009lf are rarely disturbed by discernible discrete fea¬ 
tures, especially at low latitudes, and thus this expla¬ 
nation is not a compelling one. 

L87 also introduced a condensed fraction of the 
methane within the putative cloud layer, but did not 
publish the inferred values, nor how these values were 
constrained. If the temperature profile is constrained to 
follow the wet adiabat (or the weighted average of wet 
and dry adiabats given above) then the only remaining 
variable that needs to be adjusted is the fraction of to¬ 
tal methane. There is no need to partition a fraction 
of the total into condensed form unless leaving the total 
in gaseous form would lead to supersaturation. In the 
latter case, it is reasonable to treat the excess vapor as 
condensed material. However, it is not reasonable to al¬ 
low more than a tiny fraction of the methane to be in 
condensed form. The condensed fractions reaching 10% 
or so that are noted by L87 are grossly inconsistent with 
near-IR and CCD spectral constraints because those con¬ 
straints permit only a low opacity cloud layer. 

Even a few percent of methane in condensed form 
would lead to extremely large optical depths at the cloud 
layer, which would provide very obvious spectral signa¬ 
tures that are not seen. Models of near-IR and visible 
spectra require only small optical depths, typically unity 
or less at visible wavelengths. When treated as Mie parti¬ 


cles the inferred particle size of a compact cloud near 1.2 
bars is of the order of r = 1 ;jm. The total mass per unit 
area for a given optical depth r is given by m = g prr/Q, 
where p is the density of solid methane and Q is the ex¬ 
tinction efficiency. Assuming r = 1 /zm, r = 1, Q = 1, 
and p = 0.5 g/cm 3 dCostantino and Danielsll19751 ). the 
mass density is 6.7xl0 -5 g/cm 2 . This is a factor of 
300,000 times smaller than the typical 20 g/cm 2 of total 
gaseous methane within the cloud layer. Spectral con¬ 
straints thus require that the condensed fraction must 
be so small as to play an insignificant role in the refrac¬ 
tivity profile. 

While a T(P) solution consistent with the refractiv¬ 
ity measurements is possible with no methane at all in 
the atmosphere (Model A of L87), this was rejected be¬ 
cause methane clearly plays a major role in shaping the 
spectrum of Uranus. As the mixing ratio above the 
cloud is increased, the inferred temperature must in¬ 
crease to maintain the observed refractivity, and so does 
the methane mixing ratio inferred for the deep atmo¬ 
sphere. The maximum cloud-top humidity inferred by 
L87 is 53%, which leads to a deep mixing ratio of 4%, 
although this solution did not yield the highest humid¬ 
ity level within the cloud layer. Trying to increase the 
above-cloud humidity any further results in rapidly in¬ 
creasing temperatures into the cloud layer and unaccept¬ 
able superadiabatic lapse rates. None of these profiles 
yield anything close to saturation at the altitudes where 
cloud condensation is suspected. 

2.2. Reconstructing the refractivity profile. 

With the aim of conducting a reanalysis of the occul- 
tation profile with different assumptions, we first needed 
to create a detailed refractivity profile. We began with 
the tabulation of P, T, molecular weight, number density, 
and mixing ratio published by L87. We used refractivity 
values per molecule of Kue= 0.5062, Kh 2 =0.1302, and 
Kch 4 = 1-629 (all in units of 10~ 17 cm' 3 ), which are the 
same as those used by L87 and referenced therein. We 
then used the relation 

K(z) = K H 2 fH 2 (%) + K He fHe{z) + KcH 4 fcHi ( z ) ( 2 ) 

to compute the mean molecular refractivity, where sub¬ 
scripted / values denote numeric fractions for each 
molecule. We also used the same 15/85 ratio of He to 
H 2 as L87. The total refractivity for the mixture is then 
given by 

N(z) =n(z)K(z), (3) 

where n{z) is the total number density at altitude z. Our 
computed refractivity profile versus altitude is shown in 

Fig- CD 

To validate that the refractivity profile we computed 
was consistent with the temperature profile, we inverted 
the profile as follows. We computed the pressure vs al¬ 
titude under the assumption of hydrostatic equilibrium 
using 

P(z) = P 0 +f M(z)(n(z)/N A )g(z)dz, (4) 

J Zq 

where M(z) is the molecular weight in grams per mole, 
Na is Avogadro’s number, and g{z) is the gravitational 
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Fig. 1.— Refractivity profile computed from the L87 model D 
profile. The inset provides a detailed view of the sudden slope 
change, and there horizontal lines indicate cloud boundaries for 
our D1 (dotted) and FI (dashed) structure solutions (discussed in 
Sec. 3). 

acceleration as a function of altitude, which varies from 
8.6843 m/s 2 at 1 bar (0 km) to 8.5157 m/s 2 at 240 km, 
assum ing a retrograde wind of 100 m/s (iSromovskv et al.l 
I2009D . which reduces g by only 0.23% and could be ig¬ 
nored. We started the downward integration at 240 km 
above the 1 -bar level, and took the starting pressure to 
be 2.5xl0 ~ 4 bar to match the stratospheric profile of 
L87. (The specific value of Po has an insignificant effect 
on the structure for pressures greater than a few hundred 
millibars.) Assuming an ideal gas equation of state we 
then computed T(P) from pressure and number density. 
Our T(P) profile thus constructed is compared with the 
L87 profile in Fig. 0 where we see that differences be¬ 
low the tropopause are generally smaller than 0 . 1 ° (the 
RMS deviation is 0.07 K). This provides a reasonable 
validation of our reconstructed refractivity profile, which 
we will reanalyze with different assumptions in a subse¬ 
quent section, after first validating our methane retrieval 
procedure. 

2.3. Validation of methane profile retrieval 

We next tried to reproduce the methane profile re¬ 
trievals obtained by L87. We started with the refractivity 
profile, the assumed He/H 2 ratio of 15/85, selected alti¬ 
tudes for the cloudy layer, and selected a constant rela¬ 
tive humidity above the cloud layer up to the tropopause, 
and above the tropopause assumed a constant methane 
mixing ratio equal to the tropopa use val ue, although 
later we used the 10 ” 5 upper limit of 1 Orton et all ( 1987 ). 
Within the cloud layer we forced the lapse rate to agree 
with Eq. 0 We then started at the top of the atmo¬ 
sphere and integrated downward the number density and 
pressure and iteratively solving for fcn 4 and n(z) under 



Fig. 2.— A: T(P) profile derived from refractivity profile of Fig. 
|T| (solid) using L87 tabulated results, and the L87 model D profile 
(symbols). B: the difference profile. 

the constraints of the specified methane humidity above 
the cloud, the constraints of the temperature lapse rate 
within the cloud, and the fixed mixing ratio below the 
cloud, which was taken to be the mixing ratio at the 
bottom of the cloud layer. 

Our attempt to reproduce the Model D solution is dis¬ 
played in Fig. 0 In most respects our results are barely 
distinguishable from those of L87. We obtain a deep CF 14 
VMR of 2.22% compared to their value of 2.26%, and our 
maximum CH 4 RH at the cloud base is 79.2% compared 
to their 78%. These might be brought into better agree¬ 
ment with slightly different choices for the cloud bound¬ 
aries. We did not assign any of the methane fraction to 
condensed material; we don’t know whether L87 did or 
not for this model. We take these comparisons to be ad¬ 
equate validation of our inversion technique. Note that 
the T(P) and CF 14 profiles of KT2009, which are also 
plotted in Fig. 0 have significantly higher temperatures 
and much more methane above the cloud top. 

If we increase the above cloud relative humidity as 
much as plausible, which we estimate to be about 57% 
(instead of the 53% of L87), we get results very close to 
the L87 Model F profile. We obtained a deep CH 4 VMR 
of 4.13% (instead of 4%) and a peak in-cloud humidity 
of 70% instead of 72%. The temperature profile is also 
close to the Model F profile of L87. Even at this upper 
limit, however, the above cloud methane is well below 
what was adopted by KT2009. Pushing the methane 
values above this level results in physically unacceptable 
results: the in-cloud humidity becomes lower than the 
humidity above the cloud and the sub-cloud lapse rate 
becomes more and more unstable. We would argue that 
even this upper limit leads to physically implausible pro¬ 
files because of the relatively low humidity in the cloud 
layer. Yet good fits to the spectral observations seem to 
need even more methane at these levels. In the following 
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ztop= -5.25 km, zbot= -7.75 km, rhtrop= 0.30, rhtop=0.300, exx= 1.00 fne= 0.00000 hevmr=0.1500 




ptop=1.179 bars, pbot= 1.278 bars, rhtrop= 0.30, rhtop=0.300, exx= 1.00 fne= 0.00000 He/H2=0.1765 



CH 4 (gas) mixing ratio 



Fig. 3.— A: Inverted T(P) profile for 30% RH above cloud (our Model D) compared to the L87 Model D (dotted) and the profile adopted 
by KT2009 (+). B: relative humidity profile of the inverted T(P) profile. C: Methane mixing ratio profile we inverted (solid) compared 
to the L87 Model D (dotted) and that adopted by KT2009 (+), with the saturation mixing ratio shown as the dot-dash curve, using 
our T(P) profile. D: detailed views of the methane mixing ratio profiles. The extra dashed curve is the saturation methane mixing ratio 
computed from the KT2009 T(P) profile. Note the very close agreement between our inversion and that of L87, in several cases too close 
to distinguish their difference. The horizontal dotted lines are at pressures of 1.179 bars and 1.278 bars, which correspond to altitudes of 
5.25 km and 7.75 km below the 1-bar level. 


section we show how these levels can be reached. 

3. REVISED ANALYSIS OF THE REFRACTIVITY PROFILE. 

Now that we have validated our analysis techniques, 
we apply these techniques to obtain new solutions for 
temperature and methane profiles. We first revised the 
methane mixing ratio at the tr opopause and throughout 
the stratosphere to equal the lOrton et al.l ()1987f l upper 
limit of 1x10 , and used a variable relative humidity 

between the tropopause and the cloud top, using the for¬ 
mulation 


RH(z) — RH trop + ( RH ctop — RH trop )( 1 — (5) 

Zctop)/ (ztrop Aitop)) 5 


where a tropopause humidity RH tr0 p of about 12% is 
needed to match the desired minimum methane VMR, 
and the tropopause height z c t op is taken to be 45 km. 
An exponent x=l provides linear interpolation, and a 
decrease of humidity above the cloud top that is similar 
to that adopted by KT2009. These changes made no sig¬ 
nificant difference in the plausible upper limit of methane 
mixing ratios in the vicinity of the cloud layer and in the 
deep atmosphere. We also tried adding neon to the at- 
mosphere, u sing the mixing ratio of 0.0004 suggested by 
iConrath et al.1 (119871 ). This has a very small but unde¬ 
sirable effect that reduces the upper limit on methane. 


What is really needed is a lower background molecular 
weight, which then requires increased methane to match 
the refractivity profile. The only plausible way to obtain 
a lower molecular weight is to decrease the mixing ratio 
of He, which is what we proceeded to do in the following 
fashion. Within the cloud layer we used the same formu¬ 
lation as given by Eq. [2l but decreased the He VMR (at 
all altitudes) as needed to obtain methane condensation 
within most of the putative cloud layer. Below the cloud 
layer we used the same methane mixing ratio as found 
at the bottom of the cloud layer. 

Decreasing the He mixing ratio has two beneficial ef¬ 
fects: it allows methane saturation mixing ratios to be 
attained within the layer where condensation is expected 
to occur, and it allows a higher maximum methane mix¬ 
ing ratio that is more likely to be consistent with near-IR 
and visible spectra of Uranus. Our first example solution 
(Model Dl, Fig. 0]) uses a He VMR of 0.126, and yields a 
saturated methane mixing ratio through the bottom half 
of the nominal cloud region, and a deep mixing ratio of 
2.22% for an above-cloud humidity of 38%. This is sim¬ 
ilar to Model D of L87, but is more physically plausible. 

A second example (Model FI), which attains the same 
deep mixing ratio of Model F (4%) is shown in Fig. [5] 
This solution is notable in having much more methane 
above the cloud layer than Model F, and provides a close 
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ztop= -5.25 km, zbot= -8.25 km, rhtrop= 0.20, rhtop=0.380, exx= 1.00 fne= 0.00040 hevmr=0.1260 




ptop=1.136 bars, pbot= 1.251 bars, rhtrop= 0.20, rhtop=0.380, exx= 1.00 fne= 0.00040 He/H2=0.1442 



Fig. 4.— As in Fig.[3]except that we plot our Model D1 profile, which uses a volume mixing ratio of 0.126 for He, an above cloud humidity 
of 38%, linearly interpolated to a tropopause humidity of 12%. The horizontal dotted lines are at pressures of 1.136 bars and 1.251 bars, 
which correspond to altitudes of 5.25 and 8.25 km below the 1-bar level. 


match to the CH 4 VMR profile adopted by KT2009, 
although our corresponding T(P) profile is somewhat 
cooler than their adopted profile, as needed to match 
the refractivity profile. The He VMR for this solution is 
0.1155, which is only 1.05a below the nominal value of 
0.15±0.033, given bv iConrath et al.l (I1987H . Our maxi¬ 
mum methane solution (our Model G in Table [TJ) is ob¬ 
tained for a He VMR of 0.1063, which is only 1.3a below 
the nominal value. This solution has an above-cloud hu¬ 
midity of 100%. This provides a deep CH 4 VMR of 4.88% 
and even more methane above the cloud top than that 
adopted by KT2009. 

A summary of the above model profiles and interme¬ 
diate model EF is provided in Fig. [ 6 ] and detailed pa¬ 
rameter information in Table [0 The small difference in 
tropopause temperatures (at p = 0.1 bar) is due to the 
different helium mixing ratios used in each model. The 
maximum temperatures reached near p= 2.3 bars for our 
models Dl-G are comparable to those obtained by L87 
for their models C-F. 

We now have a range of solutions that are consis¬ 
tent with occultation results, as well as being consis¬ 
tent with the expectation that methane humidity levels 
should reach saturation levels in the putative condensa¬ 
tion region, which is presumably the physical reason for 
the sudden change in refractivity slope. The remaining 
question is whether it is possible to fit the spectra of 
Uranus with a cloud particle layer in the same pressure 
regime where occultation analysis implies a cloud layer. 


We first describe our chosen spectral constraints and ra¬ 
diation transfer and fitting methods, and then proceed to 
describe the results of applying the spectral constraints. 

4. USE OF URANUS SPECTRAL CONSTRAINTS 
4.1. Spectral observations 

We chose for our spectral comparisons the spatially re¬ 
solved STIS observations made on 19 August 2002, as 
corrected and calibrated by KT2009. These spectra are 
undoubtedly the most accurately calibrated and char¬ 
acterized spectra available for the purpose. The data 
provide two spatial dimensions, with a pixel size of 0.05 
arc seconds on a side, providing 37 samples from center 
to limb, and one spectral dimension, which is sampled 
at 0.4-nm intervals, providing 1-nm resolution from 300 
nm to 1000 nm. A sample image from this cube is pro¬ 
vided in Fig. 0V These data provide a good view of the 
low latitude region sampled by the Voyager radio occul¬ 
tation experiment, although some 16 years later, which 
is a delay of 20% of a Uranus year. However, given the 
l ong rad iative time constants in the Uranus atmosphere 
(jConrath et al.llT990h , it is plausible to expect only small 
changes in cloud structure in low latitude regions. Stabil¬ 
ity at low latitudes also seems to be compatible with the 
analysis of HST images from 1994 to 2000 iKarkoschkal 
( 120011 ) . which indicate little change in albedo structure 
at these latitudes. 

In addition to spectral constraints, the KT2009 data 
provide important center-to-limb (CTL) information 
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ztop= -5.90 km, zbot= -7.75 km, rhtrop= 0.20, rhtop=0.830, exx= 1.00 fne= 0.00040 hevmr=0.1155 




ptop=1.146 bors, pbot = 1.221 bors, rhtrop= 0.20, rhtop=0.830, exx= 1.00 fne= 0.00040 He/H2=0.1306 




FeblO-125618-2011 

Fig. 5.— As in Fig. f4l except that we show our Model FI profile, which uses a volume mixing ratio of 0.1155 for He, and an above cloud 
humidity of 83%. The horizontal dotted lines are at pressures of 1.179 bars and 1.278 bars, which correspond to altitudes of 5.25 and 7.75 
km below the 1-bar level. 




Fig. 6.— A summary plot comparing the profiles of temperature (left) and methane mixing ratios (right) for models Dl, El, EF, FI, 
and G. 
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(Fig. [7)3). To reduce the effects of noise, we fit the center- 
to-limb scans to a smooth function of /i (the cosine of the 
zenith angle), as illustrated in Fig. [7)3 for a planetocen- 
tric latitude of 5° S. We then sampled that function at /i 
= 0.3, 0.4, 0.6, and 0.8. The resulting substantial noise 
reduction is readily apparent in the figure, where uncer¬ 
tainties are shown by roughly parallel dot-dash lines at 
la limits. We were able to ignore the difference between 
solar and observer zenith angles because the STIS obser¬ 
vations were taken at a very small phase angle (0.04°). 
We also were able to ignore possible longitudinal varia¬ 
tions in cloud structure because the discrete cloud fea¬ 
tures of Uranus are generally so small and of such low 
contrast that they do not obscure the CTL information 
of the zonal bands. We also know from prior experi¬ 
ence with differences between Uranus images taken at 
substantially different central meridian longitudes that 
the only substantial longitudinal I/F variations are due 
to discrete features. At the few latitudes where discrete 
features were found in the STIS data the fits easily in¬ 
terpolated across them. The selected 5° S STIS spectra 
interpolated to five cosine values are displayed in Fig. 0 
Note the strong limb darkening at short wavelengths and 
the strong limb brightening at centers of the methane ab¬ 
sorption bands at longer wavelengths. The spectrum for 
/t = 0.2, which is noisier and not well corrected because 
it is too close to the limb, was not used in our analysis. 

To further reduce noise in the observations, and to fa¬ 
cilitate model comparisons, we also smoothed the ob¬ 
served spectra to a uniform wavenumber resolution us¬ 
ing a 36 cm^ 1 boxcar. A uniform wavenumber resolution 
was chose n to fit the require ments of our Raman scatter¬ 
ing code (iSromovskvll2005al l. and the specific resolution 
is chosen to be both commensurate with the wavenum¬ 
ber shifts of the three most important Raman transitions 
and to approximately match the resolution provided by 
the STIS observations at 550 nm. The spectral smooth¬ 
ing and sampling of the fits to the CTL variations made 
noise in the observations a fairly small contributor to the 
total uncertainty as demonstrated in Sec. 4.4. 

The sensitivity of these spectra to different atmo¬ 
spheric levels on Uranus is indicated in Fig. [HI This indi¬ 
cates the penetration depth of light into an aerosol-free 
atmosphere of Uranus by the plot of pressures at which 
one-way unit vertical optical depth is reached. Pene¬ 
tration depths for individual contributions by Rayleigh 
scattering, methane absorption and hydrogen collision- 
induced absorption (H 2 CIA) are also shown. A key 
wavelength region is near 0.825 /mi, where H 2 CIA opac¬ 
ity exceeds the opacity of methane. Poor fits in this 
region are an indication of incompatible vertical distri¬ 
butions of methane and hydrogen absorptions, and may 
be a result of assuming an incorrect methane mixing ra¬ 
tio. 

4 . 2 . Radiation transfer calculations 

We used the r adiation transfer code described by 
iSromovskv! (2 005 a), which include Raman scattering and 
polarization effects on outgoing intensity. To save com¬ 
putational time we empl oyed t he accu rate polarization 
correction described bv iSromovskvl (l2005bll . After trial 
calculations to determine the effect of different quadra¬ 
ture schemes on the computed spectra, we decided to use 
14 zenith angle quadrature points per hemisphere and a 


14-order azimuthal expansion for our compact model and 
10 quadrature points for fitting models with the KT2009 
structure. To characterize methane absorption we used 
the corrected coefficients of KT2009. To model collision- 
induced absorption (CIA) of H 2 -H 2 and He-EU interac¬ 
tions we interpolated tables of absorption coefficients as 
a function of pressure and temperature that were com- 
p uted w ith a progra m provided by Alexandra Borysow 
(iBorvsow et al.ll2000l ). and available at the Atmospheres 
Node of NASA’S Planetary Data System. We assumed 
equilibrium hydrogen for most calculations but did look 
into the effects of non-equilibrium distributions, which 
are discussed in Sec. 18.21 

4 . 3 . Cloud models 

We used two distinct models of cloud structure for 
comparison purposes. The first is nearly identical to the 
model of KT2009, which we will refer to as the KT2009 
model, and the second is a modified version, which we 
will refer to as the compact cloud layer model. A plot of 
the vertical distribution of these layers can be found in 
Sec. 5.3 (Fig. 14). 

4 . 3 . 1 . The KT2009 model 

This model has four layers of aerosols, the uppermost 
being a Mie-scattering stratospheric haze layer charac¬ 
terized by an optical depth at 0.9 /im, a gamma size dis¬ 
tribution dHansenl 119711) . with a mean radius of a =0.1 
fim and a normalized variance of b =0.3. These particles 
are assumed to have a real index of 1.4, and an imaginary 
index following the KT2009 relation 

rij(A) = 0.055 exp[(350 — A)/100], (6) 

for A in nm. This haze was distributed vertically above 
the 100 mb level with a constant optical depth per 
bar. The remaining layers in the KT2009 model are 
characterized by a wavelength-independent optical depth 
per bar and a wavelength-dependent single-scattering 
albedo, given by 

za t (X) = 1 - 1/[2 + exp[(A - 290)/37]], (7) 

again for A in nm. Their adopted double Henyey- 
Greenstein phase function for the tropospheric layers 
used g\ =0.7, </2=-0.3, and a wavelength-dependent frac¬ 
tion for the first term, given by 

/i(A) = 0.94 - 0.47sin 4 [(1000 - A)/445], (8) 

which produces a backscatter that decreases with wave¬ 
length, as shown in Fig. [9] The three tropospheric lay¬ 
ers are uniformly mixed with gas molecules, with differ¬ 
ent optical depths per bar in three distinct layers: 0.1- 
1.2 bars (upper troposphere), 1.2-2 bars (middle tropo¬ 
sphere), and P>2 bars (lower troposphere). These op¬ 
tical depths are the adjustable parameters we use to fit 
this model to the observations. 

4 . 3 . 2 . The compact cloud layer model 

This model is a modification of the KT2009 model. 
The main change we made is to replace their middle 
tropospheric layer with two compact layers: an upper 
tropospheric compact cloud layer and a middle tropo¬ 
spheric compact cloud layer. This allows the possibility 
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TABLE 1 

OCCULTATION ANALYSIS PARAMETERS AND CHARACTERISTIC RESULTS FOR OUR MODELS. 


Model name: 
Tropopause RH 
Cloud top RH 
Maximum RH 
Neon VMR 
Cloud Top,km 
Cloud Bot.,km 
Cloud Top, bar 
Cloud Bot., bar 
He VMR 
AVMR/ctvmr 
H e/H 2 
Deep CH 4 VMR 
CH 4 km-am to Cld Bot. 


D 

D1 

El 

30% 

20% 

20% 

30% 

38% 

70% 

79% 

100% 

100% 

0.0 

0.0004 

0.0004 

-5.25 

-5.25 

-5.50 

-7.75 

-8.25 

-8.00 

1.179 

1.136 

1.142 

1.278 

1.251 

1.241 

0.15 

0.126 

0.122 

0.00 

-0.73 

-0.85 

0.1765 

0.1442 

0.1390 

2.22% 

2.22% 

3.24% 

0.314 

0.340 

0.576 


EF 

FI 

G 

20% 

20% 

20% 

78% 

83% 

100% 

100% 

100% 

100% 

0.0004 

0.0004 

0.0004 

-5.70 

-5.90 

-6.50 

-7.80 

-7.75 

-7.75 

1.142 

1.146 

1.153 

1.226 

1.221 

1.205 

0.1179 

0.1155 

0.1063 

-0.97 

-1.05 

-1.32 

0.1337 

0.1306 

0.1189 

3.64% 

4.00% 

4.88% 

0.610 

0.658 

0.719 




Wavelength, nm 


Fig. 7.— A: Sample image from the unsmoothed KT2009 STIS data cube at 620 nm, with grayed pixels indicating our center-to-limb 
sampling at 5° S, with grid lines at 30° intervals. B: Sample center-to-limb scans at four wavelengths (after spectral smoothing with a 
36 cm -1 boxcar), with fits shown by dot-dash lines bounding 1 a uncertainty limits. C: Sample interpolated STIS I/F spectra at 5° S 
planetocentric latitude after CTL fitting and spectral smoothing. 
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Fig. 8.— Pressure at which the vertical optical depth to space 
reached unity, shown for individual components including Raman 
scattering (dotted), Rayleigh scattering (dashed), methane absorp¬ 
tion (triple-dot-dash), and H 2 CIA (dot-dash). The pressure for all 
effects combined is shown by the dark solid curve for unit optical 
depth and by the light solid curve for r = 10. 



solid: r=0.100 fim vor= 0.30 qvols= 0.493 0.088 0.012 



Fig. 9.— Scattering properties of the stratospheric haze model 
(solid), the middle tropospheric Mie layer (dot-dash) and KT2009 
tropospheric particles (plus signs). A: extinction efficiency (Q, 
black), single-scattering albedo (vj . red) and backscatter phase 
function (P(7r), green). B: backscatter efficiency. 

of a better match to the observations if the aerosols be¬ 
tween 1 and 2 bars do not match the KT2009 assumption 
of being uniformly mixed with the gas. The upper tro¬ 
pospheric cloud (UTC) in our model is composed of Mie 
particles, which we characterized by a gamma size distri¬ 
bution with a mean particle radius of 1.2 pin and a fixed 
normalized variance of 0.1, a fixed refractive index of 1.4, 
and an imaginary index of zero. The particle radius was 
fixed at the given value because it did not vary much 


from that value in preliminary fits. For the middle tro¬ 
pospheric cloud (MTC) in our model we used particles 
with the same scattering properties as given by KT2009 
for their tropospheric particles. Both of these compact 
layers have the bottom pressure as a free (adjustable) 
parameter and a top pressure that is a fixed fraction of 
0.93 of the bottom pressure. This degree of confinement 
is approximately the same as obtained for the cloud layer 
inferred from the occultation analysis. We assumed sim¬ 
ilar confinement for the deeper layer in our model, but it 
could easily be more vertically diffuse than we assumed, 
as long as the effective pressure is similar. 

The last change we made was to replace their lower 
tropospheric layer by a compact cloud layer at 5 bars 
(the LTC), with adjustable optical depth and with the 
KT2009 tropospheric scattering properties. We found 
that this layer was needed to provide accurate fits near 
0.56 and 0.59 /im, but its pressure is not well constrained 
by the observations (the effect of varying the pressure 
can be compensated by varying the optical depth, to 
produce essentially the same fit quality). Whether this 
deep cloud is vertically diffuse or compact also cannot be 
well constrained, and thus our assumption of a compact 
cloud for this layer is a matter of convenience and is not 
compelled by observations. The wavelength dependence 
of the backscatter phase function, extinction efficiency, 
and single scattering albedo of the stratospheric haze 
are given in Fig. [9] for both the Mie particles we used 
for the putative methane layer (the UTC) and for the 
KT2009 tropospheric particles. The latter have wave¬ 
length independent optical depth, and thus the way its 
backscatter efficiency varies with wavelength is entirely 
determined by the phase function (defined by Eq. 0. 
The best-fit Mie particle size results in a smaller vari¬ 
ation in backscatter efficiency, although both decrease 
with wavelength. 

In summary, we consider two models. The diffuse one 
has the KT2009 structure, which provides a fitting stan¬ 
dard of comparison. The compact model, the main fea¬ 
ture of which is the splitting of the middle tropospheric 
layer of KT2009 into two layers, allows us to see if a 
compact layer of methane particles can provide good fits 
to the observed spectra, and to see which occultation- 
derived profile of temperature and methane mixing ra¬ 
tio provides (1) the best fit to the spectra and (2) the 
best agreement between the fit pressure for the middle 
tropospheric layer and the pressure inferred from the oc¬ 
cultation analysis. Hopefully, the best spectral fit would 
occur for the same profile that provided the best pressure 
match. As shown in the following, that is roughly what 
happened. 

4.4. Error model 

To measure fit quality we use x 2 : which requires an 
estimate of the expected difference between a model and 
the observations due to the uncertainties in both. We 
roughly characterized the uncertainty in measurements 
by the expected uncertainty in the samples of the smooth 
fits, as described previously. This uncertainty is shown 
by the purple curve in the upper panel of Fig. [101 There 
is also an overall calibration uncertainty, estimated by 
KT2009 to be 5%. However, calibration errors for spec¬ 
tra are similar to nearly wavelength-independent scale 
factor errors, which tend to cause changes in inferred op- 
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WAVELENGTH, /im 



Fig. 10.— A: Estimated contributions by measurement uncer¬ 
tainty (purple), methane coefficient uncertainty (blue), and overall 
modeling errors (red), to the combined relative I/F uncertainty for 
//=0.G (black). B: Combined Fractional error vs. wavelength at 
zenith angle cosines of 0.3 (purple), 0.6 (black), and 0.8 (green). 


tical depths of aerosols, with very little effect on their in¬ 
ferred vertical distributions, and thus we did not include 
this as an error source. KT2009 estimated relative uncer¬ 
tainties of their corrected spectra to be within 1 %, which 
were estimated by comparison of the corrected spectrally 
weighted STIS observations to band-pass filtered images. 
To this relative calibration uncertainty we added an over¬ 
all modeling uncertainty of 1 % for a combined relative 
fractional error of 1.4%, shown by the red curve in Fig. 
EH Another important source of uncertainty is due to un¬ 
certainty in the methane absorption coefficients, which 
is not easy to characterize. We assumed that the un¬ 
certainty in k had the form e(fc) = ak + eo, where we 
adopted values of a = 0.02 and eo = 5xl0 -4 km-am. 
The 2% scale-factor component (a) is roughly in agree¬ 
ment with the verification provided by the Descent Im¬ 
ager/Spectral Radio meter ( DISR) measur ements within 
Titan’s atmosphere dTomasko et al.l [20081) over the 0.5- 
0.75 pm wavelength range, though larger errors are in¬ 
dicated at longer wavelengths. The offset error is even 
less certain. KT2009 made cha nges as large as 0.01 km- 
am from previous coefficients (lKarkoschkall998t) , and the 
new values, which we use here, are likely to be much 
smaller, but by exactly what factor is unclear. 

We used a r ef ineme nt of the approach of 
ISromovskv and Frvl ([2010D in approximating the ef¬ 
fect of methane absorption uncertainty on the spectrum. 
We used a similar crude approximation, that the 
reflected intensity at any wavelength could be expressed 
as a maximum times an exponential in optical depth, 


i.e. 

I(X,p)=I 0 (p)exp(-ku/p) (9) 

where the optical depth is ku , k is the absorption coeffi¬ 
cient, and u is the absorber path amount. This directly 
leads to a fractional error in I/F given by 

SI/I=l-[I(X,p)/I 0 (p)] a+e ° /k (10) 

which differs from the model of ISromovskv and Frvl 
(|2010l) by using the maximum Iq (/i) instead of unity, and 
by including the eo /k term in the exponent. That term 
becomes dominant for small k, and expresses the fact 
that small absorptions that cause large I/F depressions 
lead to large uncertainties because the offset uncertainty 
begins to be comparable to the total absorption coeffi¬ 
cient. However, in this crude form the term becomes 
unstable when the absorption is very small and I/Iq is 
near 1, but where (because of our crude model) it is 
not close enough to 1 to make the eo /k term unimpor¬ 
tant. To counteract this instability we replace eo /k with 
eo/(A: + fcmin), where we chose /c m j n =0.01 to avoid exces¬ 
sively large errors for A < 0.5 /im. The final form of this 
model is shown as the CH 4 uncertainty in Fig. EH This 
was treated as normally distributed with the standard 
deviation as shown, which is another crude assumption. 
However, it is likely that the overall dominance of this 
error source at wavelengths exceeding about 0.6 pm is 
correct. 

Our final combined error estimate was the square root 
of the sum of the squares of the three sources. Although 
these error sources are not accurately known or well char¬ 
acterized by our noise model, they can be roughly vali¬ 
dated by comparing, for the best fit over all models and 
profiles, the minimum \ 2 value with the expected value 
for a perfect fit. Our very best fit over this spectral 
range (at this latitude) yielded approximately 339 for y 2 
instead of the value of 269 predicted by the error model. 
This would be the result of under-estimating the com¬ 
bined errors by 12 %, or by a defect in the physical model 
we are using to fit to the observations (at other latitudes 
we obtained % 2 values as low as 269). 

4.5. Fitting procedures 

To avoid the need for fitting an imaginary index in the 
methane layer (the UTC layer) we fit only the wavelength 
range from 0.55 /jm to 1.0 pm. To provide a reasonable 
computational speed while still sampling a wide range 
of penetration depths for each spectral band, we chose 
a wavenumber step of 118.86 cm ' 1 for sampling the ob¬ 
served and calculated spectrum. This yielded 69 spectral 
samples, each at four different zenith angle cosines (0.3, 
0.4, 0.6, and 0.8), for a total of 276 points of comparison. 
Our compact layer model has seven adjustable parame¬ 
ters (pressures and optical depths of the UTC and MTC 
layers, the optical depth per bar of the stratospheric 
haze, the optical depth per bar of the upper tropospheric 
haze (KT2009 referred to this as the upper tropospheric 
layer), and the optical depth of the LTC layer, leaving 
269 degrees of freedom. We fixed the LTC base pressure 
to 5 bars and the UTC mean particle radius at 1.2 pm be¬ 
cause these values were consistently obtained for a vari¬ 
ety of fits, and reducing the number of adjustable param¬ 
eters improved fit algorithm performance. Furthermore, 
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the LTC pressure is not well constrained by the observa¬ 
tions because its change can be compensated by changing 
the LTC optical depth. To fit the KT2009 model over 
the same range, we followed their approach by adjusting 
only the four dr/dP (optical depth per bar) values. We 
use a mod ified Leven berg -Marq uard t non-linear fitting 
algorithm (iSromovskv and Frvll20 1 (it) to adjust the fitted 
parameters to minimize % 2 and to estimate uncertain¬ 
ties in the fitted parameters. The uncertainty in % 2 is 
expected to be ~25, and thus fit differences within this 
range are not of significantly different quality. 

5. APPLICATION OF SPECTRAL CONSTRAINTS TO THE 

OCCULTATION SOLUTIONS 

5.1. Fit results for 5° S 

The fit results for the compact model for profiles Dl, 
El, EF, FI, and G are given in Tableland key results 
are plotted as a function of methane volume mixing ra¬ 
tio in Fig. HU The parameters of the aerosol model are 
very well constrained by the observations, with pressures 
constrained to a fraction of a percent and optical depths 
usually to within 5%. In Fig. HU we show several results 
that can be used to constrain the methane mixing ratio 
at 5° S: (A) the pressure of the upper tropospheric cloud 
(UTC) in comparison with the occultation cloud pres¬ 
sure, (B) the overall quality of the spectral fit, (C) the 
fit error at 0.825 /im, and (D), the He/H 2 ratio. 

The results for cloud pressure (Fig. HU A) show that 
the upper compact cloud (modeled as methane) is in 
best agreement with the occultation pressure range at 
a methane mixing ratio of 4.0%, but the match is still 
fairly close down to a mixing ratio near 3.5%. An even 
stronger constraint on the CH 4 mixing ratio is the fitting 
error in the region near 0.825 /mi (Fig. 11 IB j. where hy¬ 
drogen CIA exceeds methane absorption. If there is too 
much methane assumed in our model, then model cloud 
particles will need to move upward to compensate. That 
will place them relatively further above the absorption of 
hydrogen, and where the effects of hydrogen can be seen 
(as at 0.825 pm) the model will appear too bright rela¬ 
tive to the observations. Where the assumed CH 4 mixing 
ratio is too low (as for Model Dl at 5° S) the model I/F 
will appear be too low at 0.825 pm. This constraint 
clearly favors a mixing ratio of 4.5% with an uncertainty 
of roughly 0.7%. The overlap of the first two uncertainty 
ranges is 3.8-4.5%. The third constraint is overall fit 
quality (Fig. HUC). The best spectral fit (judged from 
the minimum x 2 in the middle panel) is for a methane 
mixing ratio near 3.6%, but the fit is nearly as good 
over a wide range from 3-4.9%. This is a relatively weak 
constraint that is easily compatible with the previous 
stronger constraints, which favor a mixing ratio of 4%. 
This makes our Model FI profile t he preferred pro file, 
even though it slightly exceeds the lConrath et al.1 (11987T) 
uncertainty limits of the He VMR, as shown in Fig. [TTD. 
We consider El, EF, and FI profiles plausible candidates 
for further analysis. Our preferred model (FI) leads to a 
compact methane condensation cloud very close to pres¬ 
sure level expected from occultation results and it is also 
a model that provides an excellent fit to the STIS spectral 
observations, especially in the key region where hydrogen 
CIA is significant. 

The best compact cloud model fit to STIS spectra at 
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Fig. 11.— Fit results for compact cloud layer models as a function 
of the deep volume mixing ratio of methane (points are plotted for 
models Dl, El, EF, FI, and G from left to right). A: pressures of 
the upper (presumably CH 4 ) and middle tropospheric cloud models 
compared to the cloud layer inferred from our occultation analysis 
(solid lines). B: Fit error at 0.825 ^m where hydrogen CIA is 
prominent. C: Minimum y 2 values (solid line) for fits to spectra 
at 4 zenith angles and 69 wavelengths from 0.55 to 1.0 fi m, with 
dotted lines indicating the la uncertainty range e xpected for y 2 . 
D: He volume mixing ratio (line) compared to the IConrath et all 
(119871) value (symbol with error bars). The grayed areas indicate 
rough regions of uncertainty/acceptability. 


5° S using the FI profile (shown in Fig. QUj) yielded an 
uncorrected % 2 of 342, which is significantly better than 
the best % 2 values of 409, 426, and 458 obtained by our 
fit of the KT2009 model to El, EF, and FI profiles re¬ 
spectively. Thus, we don’t have to give up anything in 
ht quality to obtain the compatibility between aerosol 
models and occultation models. In fact, the fit quality 
improves, as one would expect with two more adjustable 
parameters available. 

Our presumed methane cloud has an effective particle 
radius near 1.2 /im and a rather small optical depth, only 
about 0.3 2 at A = 0.5 um, which about half the value re¬ 
ported bv lRages et all (1 1 9911) for their low latit ude cloud, 
but close to the 0.4 value of lBaines et ~ahl (11995H for a disk 
average that weighted high latitudes most. The middle 
tropospheric cloud is considerably thicker, with an opti- 
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TABLE 2 

Best-fit parameters for compact cloud layer models at 5° S. 


OCCULTATION PROFILE: 

D1 

El 

EF 

FI 

G 

KEY PROFILE PARAMETERS: 

Cloud Top, bar 

1.136 

1.142 

1.142 

1.146 

1.153 

Cloud Bottom, bar 

1.251 

1.241 

1.226 

1.221 

1.205 

Deep Methane VMR 

2 .22% 

3.24% 

3.64% 

4.00% 

4.88% 

AEROSOL PARAMETERS: 

(dr/dP) Strat . H , bar -1 

0.329 

0.208 

0.178 

0.158 

0.136 

(d,T/dP) UJ H , bar -1 

0.01 

0.032 

0.035 

0.038 

0.031 

UjTC, bar 

1.45 

1.27 

1.24 

1.23 

1.19 

Pmtc, bar 

2.54 

1.80 

1.72 

1.68 

1.60 

P LTC , bar 

5. 

5. 

5. 

5. 

5. 

T UTC 

0.444 

0.328 

0.324 

0.322 

0.324 

T MTC 

1.33 

1.16 

1.23 

1.28 

1.42 

T LTC 

0.01 

2.30 

2.92 

3.45 

4.75 

r Strat.H , M m 

0.1 

0.1 

0.1 

0.1 

0.1 

r UTC, Atm 

1.2 

1.2 

1.2 

1.2 

1.2 

FIT QUALITY: 

X 2 (total) 

435.6 

340.5 

338.5 

342.3 

351.9 

x 2 (0.825 /an) 

47.7 

11.0 

5.1 

2.2 

0.6 

Fit error at 0.825 fim, fi= 0.8 

-4.54 

-2.14 

-1.34 

-0.71 

+0.54 


Note: The uncertainty in x 2 is ~25 and thus fits differing by less than this are not of significantly 
different quality. The stratospheric haze base extends upward from 100 mb, the upper tropospheric 
haze (UTH) extends from 900 mb to 100 mb, The upper tropospheric cloud (UTC) and middle 
tropospheric cloud (MTC) extend from the bottom pressures, tabulated here, to 0.93 times those 
pressures. The lower tropospheric cloud is bounded by the tabulated base pressure and 0.98 times 
that pressure. For these fits, the particle radii of the Mie layers were held fixed, as was the pressure 
of the lower tropospheric cloud. The fit error at 0.825 /im is the ratio of (model-measured) to 
estimated uncertainty. 

cal depth of 1.23 for the optimum methane profile. The 
much deeper lower tropospheric cloud is generally even 
thicker, except for the model D1 profile for which this 
cloud is not even needed. 

5.2. Layer contributions 

The relative roles played by our five model layers in 
creating the observed spectral characteristics are illus¬ 
trated in Fig. UTTll-F. which displays the effect of remov¬ 
ing each layer from the model spectrum. The distinctly 
different contributions of the various layers is what allows 
the model parameters to be so well constrained by the 
observations. We see in B that the stratospheric haze 
(layer 1) serves to reduce the I/F for A < 0.6 (im and 
provide a slight (5-10%) boost to the I/F in the center 
of the deep absorption bands at longer wavelengths. As 
shown in C, the upper tropospheric haze layer (layer 2) 
makes a similar contribution but without as much short¬ 
wave absorption. The methane cloud (layer 3, panel D) 
contributes almost nothing at wavelengths less than 0.55 
fim but provides a significant contribution to minima in 
the intermediate absorption bands and especially to the 
shoulders of the strong absorption bands. The effect of 
the middle tropospheric cloud (layer 4, panel E) is seen 
mainly in the strong contributions to the I/F peaks, and 
also in contributing a small absorption peaking near 0.4 
fim. The lower tropospheric cloud (layer 5, panel F) con- 
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Fig. 12.— A: STIS 5° S spectra (solid lines) at four zenith an¬ 
gles, sampled at a wavenumber spacing of 118.86 cm -1 , compared 
to model calculations (points) for the model that best fits the spec¬ 
trum when the FI profile of temperature and methane mixing ratio 
is assumed. B: Ratio of Model to observed I/F values. C: Differ¬ 
ence of model and observed values divided by the expected uncer¬ 
tainty at each wavelength sample. Note that the larger fractional 
errors that occur for A >0.78 //m are within the expected range of 
uncertainties. 
















14 


Sromovsky et al. 


tributes mainly to the longer wave window regions, and 
is unique in providing key contributions at 0.55 and 0.59 
/im. Note that, although the model fit was only con¬ 
strained by measurements from 0.55 to 1.0 /im, it pro¬ 
vides a relatively good fit to the shorter wavelengths as 
well (panel A). The discrepancies between model fit and 
measurements between 0.3 and 0.4 /im are only about 
10 % and could be largely removed by slight increases in 
the lower tropospheric cloud absorption and slight de¬ 
creases in the stratospheric haze absorption in this re¬ 
gion. It might also be possible to introduce some of the 
extra absorption needed at small zenith angles i nto th e 
methane layer itself, as suggested bv lRaees et al.l (jl991i h 

5.3. Comparison of vertical structures 

The vertical structure of our compact cloud model 
and the diffuse vertical structure of the KT2009 model 
are compared directly in Fig. 1141 which displays opti¬ 
cal depth per unit pressure versus pressure (A) and cu¬ 
mulative optical depth versus pressure (B). There is a 
very large difference in the vertical distribution of cloud 
particles, with our compact model providing high opac¬ 
ity concentrations in narrow pressure ranges. But when 
compared on the basis of cumulative opacity we see that 
the KT2009 model looks like a vertically averaged ver¬ 
sion of our results. The compact layer with the strongest 
effect on the CCD spectrum is our layer 4 (MTC), as 
shown in Fig. [T3] The dashed line in Fig. ITTlV is f or op¬ 
tical d epth per bar at 1.6 p m from Fig. 5 of llrwin et al.l 
( 12010 H with pressures scaled by the factor 1 . 8 / 2.6 to ac¬ 
count for the decrease indicated in their Fig. 8 , when 
they used a vertical methane profile similar to the L87 
Model F. Their relative vertical distribution of opacity is 
crudely consistent with that of the other models, when 
vertical smoothing is taken into account. However, the 
quantitative values of their opacities per unit pressure 
are hard to compare with our (or KT2009) results be¬ 
cause of very different assumptions made about particle 
scattering properties, including the extreme assumption 
that particle scattering properties were independent of 
latitude and altitude. 

While the top compact layer in our model is consis¬ 
tent with condensation of methane, the composition of 
the deeper layers is highly uncertain. Possible parent 
gases are H 2 S and NH 3 , but condensation layers at 1.7-3 
bars would imply strong depletions from their expected 
deep mixing ratios, which is plausibly an effect of the 
formation of a much deeper cl o ud o f NH 4 S H particles, 
as dis cussed by Ide Pater et al.l (11989(1 and iFeelev et al.1 
(fl99l . 

6 . LATITUDE DEPENDENCE OF METHANE ON URANUS 

We first show that the occultation solution that is most 
consistent with low latitude spectra is not consistent with 
high-latitude spectra and then discuss how that profile 
can be modified in physically reasonable ways and which 
modifications provide the best fits at high latitudes. 

6 . 1 . Fit results assuming latitude-independent methane 

EF and FI model temperature and methane profiles 
yielded high quality spectral fits not only at the occulta¬ 
tion latitude, but also for latitudes from 30° S to 20° N, 
as indicated in Fig.lTSlV for the FI fits. Over this latitude 


range the fit quality was generally very high, reaching y 2 
values close to those expected from the error model, and 
about half were notably better than we obtained at 5° 
S. However, between 30° S and 45° S fit quality began 
to deteriorate dramatically, with y 2 growing to ~1100 
by 60° S, which is four times the expected value. At the 
same time, at 0.825 pm, where hydrogen CIA is an im¬ 
portant contributor, the y 2 contribution (Fig. HsBl and 
the signed fit error (Fig. [I5p) remained relatively flat 
over the 30° S to 20° N region, but grew dramatically 
from 30° S to 60° S, with the signed error reaching 7 
times the expected uncertainty and the y 2 contribution 
reaching more than 40 times the expected value. Thus, 
there is little question that the methane mixing ratio at 
high latitudes is quite different from what it is at middle 
latitudes, a result already established by KT2009. To get 
a more realistic picture of the high latitudes, we need to 
estimate what the actual methane mixing ratio profile is 
at these latitudes. In the next section we discuss plausi¬ 
ble ways to modify our baseline (FI) profile so that we 
can investigate this issue. 

6 . 2 . Construction of depleted methane profiles 

The occultation places no direct constraints on the 
temperature or methane profiles at other latitudes. How¬ 
ever, there are some reasonable physical constraints to 
consider. First, the lack of significant latitudinal gradi¬ 
ents in the temperatures inferred from Voyager 2 IRIS 
observations for P>~150 mb dHanel et al.l 11986 1 sug¬ 
gests that the occultation thermal profile is a good ap¬ 
proximation at all latitudes on Uranus. A second con¬ 
straint has to do with an extension of the relationship 
between the vertical wind shear and horizontal temper¬ 
ature gradients. The variation of mixing ratio with lati¬ 
tude, assuming the same pressure-temperature structure 
at all latitudes (as assumed by KT2009 as well), leads 
to a variation in density on constant pressure surfaces, 
which is similar in effect to horizontal temperature gradi¬ 
ents^ Both kinds of gradients lead to vertical wind shear 
( Sun et al.lll99lH . a consequence of geostrophic and hy¬ 
drostatic balance. If such latitudinal variations occurred 
through great atmospheric depths, this would lead to 
great differences in cloud level winds and conflict with 
the observed wind structure of Uranus. Thus we expect 
that where methane is depleted, it is not depleted over 
great depths. And, as indicated by KT2009, the spec¬ 
tral observations do not require that methane depletions 
extend to great depths. In fact, we will later show that 
the spectral constraints favor relatively shallow methane 
depletions. 

For trial purposes we created several types of modi¬ 
fications of the FI profile. We first defined an upper 
tropospheric mixing ratio and a set of transition pres¬ 
sures, Pi and P 2 . for P > P 2 we kept the mixing ratio at 
the FI deep value of 4%. For P < Pi we set the mixing 
ratio to the upper level value, except where that value 
exceeded the FI value above the methane condensation 
level, at which point we reverted to the FI model pro¬ 
file. Between the two transition pressures we interpolated 
between the upper troposphere and deep values. Sam¬ 
ple depleted profiles of this type are shown in Fig. [TUK. 
These profiles are consistent with methane condensation 
at somewhat lower pressures than for the FI profile. A 
depleted profile of this type might be created by descent 
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Fig. 13.— A: Comparison of measured STIS spectrum (red) at 5° S with a compact model fit (black), with spectra shown at a solar 
zenith angle cosine of 0.6 and ratios shown for all four cosine values); B-F: Comparisons of the model spectrum with spectra for same 
model with layers 1-5 removed respectively. Layer 3 (D) is the methane cloud layer. Although these spectral sensitivities were calculated 
for the EF profile, they also closely represent results for El and FI profiles. 


of the low-mixing-ratio gas from slightly above the 1.2 
bar level downward to higher pressures, but to limited 
depths. 

We also considered “proportionally descended gas” 
profiles in which the model FI mixing ratio profile a(P) 
was dropped down to increased pressure levels P'(a) us¬ 
ing the equation 

P' = P x[l + (a/a d ) vx (P d /P cb - 1)] (11) 

for P tr < P < Pd, 

where Pd is the pressure depth at which the revised mix¬ 
ing ratio a! = a(P') equals the uniform deep mixing ra¬ 
tio ad , P c b is the cloud bottom pressure (which is where 
the unperturbed profile departs from the uniform deep 
value), P tr is the tropopause pressure (100 mb), and the 
exponent vx controls the shape of the profile between 100 
mb and Pd- A sampling of profiles of this type are shown 
in Fig. I16B . The pr ofiles with vx = 1 are sim ilar in form 
to those adopted bv lKarkoschka and Tomaskol 0201 ID . A 


profile of this type is consistent with the descending gas 
beginning at lower pressures than for our initial set of 
models. In both cases we would expect methane conden¬ 
sation to be inhibited by the downward mixing of up¬ 
per tropospheric gas of low methane mixing ratios. How 
these various types of profiles affected fits to the 45° S 
and 60° S spectra is discussed in the next section. 

6.3. Constraining methane profiles at 6CP S and f5° S 

We used our 5-layer model, allowed adjustment of all 9 
parameters, and tried CH 4 mixing ratios as low as 0.3% 
in the depleted region, and transition depths ranging 
from 1.5-2 bars to as deep as 9-10 bars. The bottom of 
the transition region is where the FI deep mixing ratio 
of 4% is reached. The results were evaluated by com¬ 
paring their overall fit quality (how low their y 2 values 
were) and how big their fitting errors were at the key 
wavelength of 0.825 p m, where H 2 CIA is an important 
contributor. The results are summarized in Fig. [lT] for 
depletion depths of 1.5-2 bars (only for 45° S), 2-3 bars, 
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Fig. 14.— Optical depth per bar (A) and cumulative optical depth (B) for a vertically diffuse KT2009 model (dotted) and out compact 
layer model (solid) fit to 5° S STIS spectra using the FI vertical profile of temperature and methane mixing ratio. For the KT2009 model 
we fit the four optical depths per bar, and used the KT2009 values for pressures and particle scattering properties. For our model we fit 
two optical depths per bar for the upper two layers and the pressures and optical depths of the three compact layers, except for the fixed 
5 bar pressure we used for the bottom layer. The dashed line in A is optical depth per bar at 1.6 fim from Fig. 5 of I Irwin et al.l <120101) 
with pressures scaled by the factor 1.8/2.6 to account for the decrease indicated in their Fig. 8, when they used a vertical methane profile 
similar to the L87 Model F. 



-60 -40 -20 0 20 

Planetocentric Latitude, ° 


Fig. 15.— Spectral fit quality measured by x 2 over the entire 
range of the fit (A) and x 2 f° r the 4 view angles at the 0.825 /^m 
H 2 CIA wavelength (B), and fit error at 825 nm and /z=0.8 (C), all 
as a function of planetocentric latitude, using the FI temperature 
and methane profiles. Filled circles in A denote results of fitting 
the KT2009 structure model. The horizontal dotted lines in A and 
C indicate expected (central) and la uncertainty bounds. 


4-6 bars, and 9-10 bars. At neither latitude is it possible 
to obtain a good fit with deep depletions that extend to 
9-10 bars. Although a low hydrogen fitting error can be 
obtained from deep depletions, the overall fit quality is 
much better for shallow and very strong depletions. At 
both latitudes lower mixing ratios lead to better fits, as 
long as the depletion depth is reduced as well. The mini¬ 
mum hydrogen error moves to lower methane mixing ra¬ 
tios as the depletion depth is decreased and the depletion 
factor is increased. It is apparent that the best combi¬ 
nation of overall and CIA fits requires very low methane 
mixing ratios, of the order of 1% or even less, and shal¬ 
low depletion depths, down to only 1.5-2 bars for 45° S, 
and perhaps 2-3 bars at 60° S, which are illustrated in 
Fig. ITfiR. These perturbations of the base profile are of 
such a small depth that the vertical wind shear caused 
by the resulting horizontal density gradients would not 
likely produce significant perturbations of the zonal wind 
field, though exact calculations would be needed to verify 
that. The depletion of methane at high latitude helps to 
explain the differe nce between prio r estimates of 4% b y 
iRaees et al.1 (119911) and 1.1-2.3% bv lBaines et al.1 (119951) . 
the latter weighting high latitudes more than the former. 

A limited fitting exploration using more physically ap¬ 
pealing methane profiles of the type shown in Fig. I16B 
did not yield significantly better fits than those shown in 
Fig. [HIV. We did find constraints on the parameters of 
such profiles, however, with vx = 1 yielding better fits 
than either vx = 0.5 or vx = 2, and P p d = 2 — 3 bars 
preferred at 45° S and P p d =4 — 6 bars preferred at 60° 
S. 

These results present a picture in which dry upper 
tropospheric gas descends at high latitudes after being 
dried out by condensation of methane in upwelling gas 
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Fig. 16. — Sample profiles of depleted methane used to fit high- 
latitude STIS spectra. The undepleted base profile labeled by its 
valid latitude range of 30° S to 20° N is from Model FI. The 
profiles in A provided the best fits for the indicated latitudes. None 
of the B profiles yielded any better fits (see text). 

at low latitudes, with upwelling and descending regions 
connected by meridional cells as illustrated in Fig. [18] 
The general nature of this circulation is similar to that 
described by KT2009 except that they don’t provide for 
the existence of a methane condensation cloud, which 
seems to us a necessary mechanism to dry out the as¬ 
cending gas. The high latitude descending gas appears 
to be very dry and to descend only a small distance below 
the original methane condensation level. However, this 
descent should inhibit the formation of methane clouds 
at these latitudes, and if any do form they would have 
to form at much lower pressures. When we use five lay¬ 
ers, our layer 3, which used to be found at the methane 


condensation level, is now found much deeper in the at¬ 
mosphere and thus must be composed of some other sub¬ 
stance, perhaps part of a somewhat extended H 2 S cloud 
layer. 

7. LATITUDE DEPENDENCE OF CLOUD STRUCTURE 

Our cloud pressure and optical depth parameters as a 
function of latitude are given in in Table [3] and key re¬ 
sults are plotted in Fig. [12] The results for latitudes from 
30° S to 20° N were all obtained using our FI profile of 
methane and temperature. For 60° S we used a methane 
profile depleted to 0.9% with a transition to 4% between 
3 and 4 bars. At 45° S, we used a depletion to 0.6% with 
a transition to 4% between 1.5 and 2 bars. All three of 
these methane profiles are shown in Fig. I16A . Looking 
at the low-latitude region (30° S to 20° N), we see that 
the upper compact layer (the putative methane layer, 
also known as the UTC or layer 3) changes pressure only 
slightly, remaining consistent with methane condensation 
at the position of that layer, while its optical depth in¬ 
creases slightly towards the equator before beginning a 
significant decline into the northern hemisphere. Some¬ 
what more modulation is seen in the pressure of the 
dominant middle tropospheric cloud (MTC, layer 4) and 
there is a substantial declining trend in optical depth to¬ 
wards the north, though the rate of decline decreases near 
the equator. Though these declines may seem modest in 
the logarithmic plot, they are 7-10 times the typical 5% 
fitting uncertainty. The upper tropospheric haze reaches 
a peak opacity near the equator, as noted by KT2009. 

The prominent bright band near 45° S that is observed 
in Uranus images at wavelengths of intermediate absorp¬ 
tion (0.1-0.4 /km-arn, as shown in Fig. 23 of KT2009) is 
also pro minent in near-IR H- and J -band images, as in 
Fig. 7 of iSromovskv and Frvi (2007). According to Fig. 
d this band is in a region where cloud pressures have 
increased in both of the upper compact layers while the 
upper tropospheric methane mixing ratio has decreased. 
It is clear that the upper compact layer at 45° S and 
60° S cannot be the putative methane ice layer, because 
the pressure levels are much too high to permit conden¬ 
sation. It is conceivable that in this region the upper 
tropospheric cloud is actually what we called the mid¬ 
dle tropospheric cloud at low latitudes. And the much 
deeper second compact layer may be a new cloud layer 
formed by the upwelling produced by the lower circu¬ 
lation branch in Fig. [TS] This would suggest that the 
boundary between the upper and lower branches may be 
in the 1.5-1.7 bar region. 

The bright band contrast relative to 60° S is a result 
of the deeper cloud being much deeper at 60° S and hav¬ 
ing a much lower optical depth. The contrast relative 
to lower latitudes results from a combination of reduced 
methane opacity as well as the extra cloud opacity in the 
2 -bar region, which apparently more than compensates 
for the absence of the 1.2-bar cloud layer and the signif¬ 
icantly reduced opacity at 1.5 bars. At 60° S the middle 
tropos pheric cloud reaches about the same pressure as 
I Baines et al.l ( 1995 ) inferred for the semi-infinite cloud 
in their model, which is heavily weighted towards high 
latitudes. 


8 . DISCUSSION 
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Fig. 17.— Overall fit quality (A) and fit quality at 0.825 /rm (B) versus depleted mixing ratio for planetocentric latitudes 60° S (Left) 
and 45° S (Right). The transition pressures between upper depleted regions and the deep region at 4% are noted in the legend. 


TABLE 3 

Fit parameters at 10 latitudes, using Model FI profiles for all but 45° S and 60° S. 


Centric 

Latitude 

Pure 

bars 

Pmtc 

bars 

(. dr/dP) 

Strat.Haze 

(. dr/dP ) 

UT. Haze 

T UTC 

T MTC 

T LTC 

x 2 

20° N 

1.185 

1.565 

0.230 

0.0000 

0.163 

0.86 

2.74 

372 

10° N 

1.175 

1.595 

0.334 

0.0019 

0.213 

1.17 

4.16 

269 

4° N 

1.180 

1.558 

0.148 

0.0363 

0.269 

1.15 

3.83 

394 

0 ° 

1.192 

1.585 

0.185 

0.0531 

0.291 

1.22 

3.89 

351 

5° S 

1.225 

1.676 

0.158 

0.0375 

0.322 

1.28 

3.45 

342 

10° S 

1.224 

1.649 

0.280 

0.0098 

0.297 

1.40 

4.68 

281 

20° S 

1.222 

1.622 

0.262 

0.0043 

0.248 

1.48 

5.48 

282 

30° S 

1.192 

1.500 

0.075 

0.0246 

0.265 

1.53 

5.85 

311 

45° S 

1.533 

2.074 

0.234 

0.0181 

0.393 

2.24 

21.4 

435 

60° S 

1.589 

3.131 

0.206 

0.0493 

0.445 

0.67 

23.8 

444 


Note: Fixed parameters (except at 45° S and 60° S) included the lower cloud base 
pressure (5 bars), the particle radius of the upper tropospheric cloud (1.2 /im ). At 
45° S and 60° S, we used vjjtc = 1-1 and 1.34 fim , and Pltc = 8 and 11 bars, 
respectively. 


8 .1. Occultation failure to detect the middle 
tropospheric cloud 

Because the occultation probes to higher pressures 
than the pressures we infer for the dominant middle tro¬ 
pospheric cloud, one might wonder why the occultation 
did not detect this layer if it is indeed compact and cre¬ 
ated by condensation? That turns out not to be a con¬ 
tradiction. If the mixing ratio of the condensible is low 
enough, then the change in molecular weight might actu¬ 
ally be very small, and the resulting refractivity change 


might be too small to be detected. For example, for H 2 S 
to form a cloud base near 1.7 bars its volume mixing ratio 
would need to be ~10 - ' (^0.003 times the solar mixing 
ratio). The maximum change in molecular weight in¬ 
duced by condensation of that much H 2 S would then be 
about 3xl0 -6 , which would be completely undetectable 
in the occultation measurements, even if it occurred in 
a thin layer. But could such a cloud actually produce 
enough condensible material to produce the needed re¬ 
flectivity of the cloud? The answer is yes. If half the 
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Fig. 18. — Schematic of methane depletion through condensation 
at low latitudes as a drying process, poleward transport, then de¬ 
scent of dry gas, finally mixing with moist gas on the return flow. 
How far beyond 60° S the depletion extends towards the pole, or 
whether northern high latitudes are depleted is unknown. Also 
unknown is how deep the return flow extends and whether or not 
there is a deep meridional flow in the opposite sense. 
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Fig. 19.— Cloud model parameters as a function of latitude, from 
Table [3] The Upper Tropospheric cloud from 30° S to 20° N is 
consistent with a methane cloud. There appears to be no methane 
cloud at 45° S and 60° S. The largest effect on the STIS spectra is 
produced by the Middle Tropospheric Cloud. Th e filled cir cle in A 
marks the top of the optically thick cloud of IBaines et all (1995), 
plotted at the sub-earth latitude of the disk-integrated observations 
they analyzed. 


H 2 S gas in a 2-km thick layer were to condense into 1- 
/jm droplets, we estimate that the resulting cloud would 
have about 1 optical depth. Convective processing of 
additional gas could increase its opacity to much higher 
values. NH 3 is a much less suitable candidate because 
its corresponding mixing ratio would be three orders of 
magnitude smaller as would the mass available for cloud 
particle formation. 

8 .2. Effects of para-hydrogen variations 

All our calculations described so far have assumed equi¬ 
librium hydrogen, for which the para and ortho fractions 
are in local thermodynamic equilibrium. Since these two 
forms of hydrogen have different spectral features, this 
assumption deserves some consideration. The equilib¬ 
rium para faction for our FI profile is displayed in Fig. 
1201 along with several non-equilibrium curves, which we 
discuss in the following. In the deep atmosphere high 
temperatures produce a para fraction of 0.25 (referred to 
as normal hydrogen), but above the methane condensa¬ 
tion level the para fraction increases dramatically. By 
vertically mixing deep atmospheric gas up to the cloud 
level and higher, a dramatic change in the para fraction 
can be produced, if done on a time sca le short compared 
to the relaxation time. 1 Baines et al.l (119951) considered 
mixtures of equilibrium and normal hydrogen, as might 
be produced by vertical mixing of the high temperature 
para distribution with the locally equilibrated distribu¬ 
tion. Their results were consistent with mixtures that 
were no less than 85% equilibrium (and thus 15% nor¬ 
mal). KT2009 found that the difference between 100% 
and 85% equilibrium had a smaller effect than their un¬ 
certainties. They also stated that 50% equilibrium hy¬ 
drogen would shift their results toward lower aerosol 
opacities and lower mixing ratios, but not change their 
results about latitudinal variations unless the the equi¬ 
librium fraction itself varied with latitude, which they 
argued against on the basis of smooth latitudinal varia¬ 
tions in the observed hydrogen absorption regions. 

Mixtures of equilibrium and normal hydrogen cannot 
produce the very different para fraction profiles that can 
be produced by the descending or rising branches of our 
putative methane circulation. For convenience, we pa¬ 
rameterize the rising or descending gas by Pfac, which 
is the factor by which pressure at the midpoint of the 
equilibrium profile is displaced. In Fig. I2UR. we show 
para fractions that result when gas is vertically trans¬ 
ported upward (Pfac=0.5) or downward (Pfac=1.5 and 
Pfac=3.0), assuming that the transport time is much less 
than the relaxation time to reach local equilibrium. Note 
that mixtures of equilibrium and normal hydrogen always 
decrease the para fraction, while vertical gas transport 
can dramatically increase or decrease the local para frac¬ 
tion. How changes in the para profile affect the Uranus 
spectrum can be seen in Fig. [5T] for two profiles of the 
type shown in Fig.EDlA.: a descending gas case (Pfac=2.5) 
and an ascending gas case (Pfac=0.5). The effect is lo¬ 
cal, mainly between 0.77 and 0.86 /im, and is roughly 
the same percentage over a wide range of view angles, 
reaching a maximum difference of about 15%. For the 
case shown in Fig. UTill. where there is descending gas 
above 1.5 bars, and ascending gas below (deeper than) 
1.5 bars, which is more consistent with high-latitude 
motions shown in Fig. 1181 the spectral difference from 
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Fig. 20.— A: Profiles of para hydrogen assuming local thermodynamic equilibrium (solid), mixtures of equilibrium and normal hydrogen 
(dotted), and vertically displaced hydrogen before new equilibrium conditions are established (dot-dash). B: Example of a mixed rising- 
descending profile (dot-dash) compared to the equilibrium profile. 



Fig. 21.— A: Model spectra computed for ascending (red) and 
descending (black) gas profiles at view angle cosine =0.6. B: Ratio 
of the two spectra at four different view angle cosines. 


the equilibrium case is only a few percent, and not de¬ 
tectable within current uncertainties. The fact that we 
did not find evidence in the spectra for significantly al¬ 
tered para profiles is consistent with either such a com¬ 
bined upwelling/down-welling profile or no significant de¬ 
parture from equilibrium hydrogen. However, this issue 
may be worth further investigation, especially if methane 
absorption coefficients become better known. Part of 
the problem making use of this constraint is that we are 
working with a part of the spectrum where methane ab¬ 
sorption is weak and thus currently is significantly un¬ 
certain. 

8 . 3 . Plausibility of a reduced He VMR on Uranus 

The Voyager Infrared Radiometer Interferometer 
and Spectrometer (IRIS) and Radio Science Sub¬ 
system (RSS) together yielded a determination of 
Jupiter’s He/H 2 ratio equal to 0.114±0.0 25, w hich 
is a weighted mean of results by iGautier et al.1 


(119811) and IConrath and Gieraschl (119841) . However, 
the in situ measurements of the Galileo Probe con- 
tradicted this value, instead obtaini ng 0.157T0.004 
(Ivon Zah n et al.l ll998i) and 0 .156T0.006 (jNiemann et al.1 
ll998f )7 lConrath and Gautieil (|2000D reviewed possible ex¬ 
planations for this discrepancy and ruled out IRIS cal¬ 
ibration errors, but were unable to identify any plau¬ 
sible error in the occultation measurements or analy¬ 
sis, although a detailed error analysis of the occulta¬ 
tion measurements was not conducted. This discrep¬ 
ancy also cast suspicions on the very low He/H 2 ratio 
obtained by IRIS-RSS analysis for Saturn, and led to an 
upward revision of the Saturn ratio from 0.034±0.024 to 
0.11- 0.16, based o n a rather uncertain IRIS-only analysis 
(jConrath and GautieilI2000i ). This also raises questions 
about the occultation results for Uranus, suggesting per¬ 
haps that the IRIS-RSS determined ratio is actually too 
low for Uranus, rather than our suggestion that the ra¬ 
tio is too high. But an even higher ratio, with the same 
refractivity profile would lead to very serious amplifica¬ 
tion of inconsistencies with near-IR and CCD spectra 
of Uranus that have already been noted. Furthermore, 
because the specific error that caused prior discrepan¬ 
cies is not known, it is not possible to predict with any 
confidence how it might affect the Uranus analysis, or 
even what direction it might take. The change in re¬ 
fractivity per molecule that results from changing the 
He VMR from 0.15 to 0.11, as computed from Eq. [2] is 
from 4.496xl0~ 18 /cm 3 to 4.647xl0 -18 /cm 3 , which is a 
change of 3.3%. We do not know if an error of this mag¬ 
nitude is even remotely conceivable for the occultation 
measurements. 

8 . 4 . The reasons for higher cloud pressures obtained 
from Near-IR observations 

The optical depth of the methane cloud is relatively 
small and the influence of the middle tropospheric cloud 
is much more significant because of its greater opacity. 
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Fig. 22.— Comparison of methane profiles from our EF and FI 
models with those used by KT2009, L87, and llrwin et all 11201 11 1, 
the latter following the L87 profile above the point of intersection 
of their assumed deep 1.6% VMR with that profile. 

That and the limited vertical resolution inherent in all 
the spectral observations, will lead inversion methods 
such as the NEMESIS routine used in by Irwin and col¬ 
leagues (|Irwin et al . 2010, 20fH) to find peak opacity con¬ 
tributions at pressures deeper than that of the methane 
cloud. In addition, these and most prior modelers of 
near-IR spectra assumed profiles of methane with much 
smaller column amounts of methane, often by a factor of 
two or more (Fig. l22l) . which leads to much higher pres¬ 
sures for the level of significant aerosol scattering. The 
dependence of derived pr essure on assu med methane was 
clearly demonst rated bv iSromovskv and Frvl (l2007f) and 
by llrwin et al.1 (l2010l) . As noted previously, when Irwin 
et al. switched from Model D to the Model F profile, 
the pressure of their peak opacity concentration moved 
from 2.5 to 1.7 bars. Since our Model FI has more col¬ 
umn abundance of methane down to the 1.2 bar level 
than even the L87 Model F, we expect similar pressures 
would be obtained using our recommended profile. 

9. SUMMARY AND CONCLUSIONS 

After reanalysis of the L87 radio occultation profiles of 
refractivity vs altitude, and fitting compact cloud layer 
models to STIS spectra as corrected and calibrated by 
KT2009, we reached the following main conclusions: 

1. By decreasing the stratospheric He mixing ratio 
from its nominal value of 0.15 by 1-1.3 times its 
uncertainty it is possible to achieve methane satu¬ 
ration within the layers suspected to have conden¬ 
sation and to achieve increased methane humidi¬ 
ties above the condensation level, even exceeding 
the high values adopted by KT2009 at low lati¬ 
tudes, which are inconsistent with the original L87 
profiles. The maximum deep mixing ratio that we 
could obtain within reasonable physical constraints 
was 4.88% (for our model G). 

2. A five-layer cloud model in which the bottom two 


diffuse layers of the KT2009 model are split into 
three compact layers, when constrained by STIS 
spectra at 5° S, yield best-fit pressures for the top 
compact layer in excellent agreement with the loca¬ 
tion of the occultation cloud layer for profile mod¬ 
els with deep CH 4 mixing ratios between 3.2 and 
4.5%, with the best compromise fit being obtained 
at 4% (for our Model FI), although this fit is not 
significantly better than for the other models in 
this range. 

3. As judged by fitting errors at 0.825 //m, where H 2 
CIA exceeds methane absorption, the best compat¬ 
ibility between methane and H 2 CIA is obtained for 
a methane mixing ratio of 4.5% with an uncertainty 
range of about 0.7%. 

4. When we consider constraints of upper cloud pres¬ 
sure, fit error at 0.825 /im , overall fit quality, 
and t he uncertainty limits of the iConrath et al.1 
(ll987f l helium abundance, our best compromise es¬ 
timate for the deep methane mixing ratio at 5° S is 
4.0±0.5%, and the corresponding preferred vertical 
temperature and methane profiles are embodied in 
our FI Model. 

5. Our five-layer cloud model, using the EF and/or FI 
profiles, can also provide excellent fits between lat¬ 
itudes of 30° S and 20° N, with relatively latitude- 
independent cloud pressure boundaries, but gener¬ 
ally decreasing optical depths from south to north. 
The main exception is the upper tropospheric haze, 
which shows a strong peak near the equator, as 
noted by KT2009. The lower tropospheric cloud is 
found to have the greatest opacity, but the 1-2 op¬ 
tical depth (at 0.5 /im) middle tropospheric cloud 
(in 1.5-1.7 bar range) has the most significant effect 
on the observed CCD spectrum. The pressure of 
the 0.15-0.3 optical depth methane cloud layer did 
not vary much with latitude, staying within 2 % of 
its 1.2 bar mean. 

6 . Trying to fit our model at high southern latitudes 
made it clear that the methane mixing ratio profile 
could not remain the same as we used in the low lat¬ 
itude observations. Not only did overall fit quality 
deteriorate, but the fit quality at 0.825 /im wave¬ 
length of H 2 CIA got suddenly very bad, and the 
direction and size of the error at that wavelength 
indicated a significant lowering of the methane mix¬ 
ing ratio, as first pointed out by KT2009. 

7. We created depleted methane profiles in which de¬ 
pletions reached a limited depth at which point the 
mixing ratio transitioned to the same deep value 
of 4% used at low latitudes. We did not expect 
a depletion at great depths because the horizontal 
density gradients that would be generated would 
induce problematic vertical wind shears. Overall 
and 0.825 /im fit quality was found to be mini¬ 
mized with mixing ratios less than 1 % and deple¬ 
tion depths of 3-4 bars at 60° S and only 1.5-2 bars 
at 45 ° S. 

8 . Creating the observed methane depletion is plausi¬ 
bly the result of upwelling at low latitudes, where 
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the Uranus atmosphere is dried out by condensa¬ 
tion of cloud particles, and subsequent poleward 
transport to and descent of the dried-out atmo¬ 
sphere at high latitudes down to levels of a few 
bars, with greater descent nearer to the pole. This 
would inhibit methane condensation at high lati¬ 
tudes and perhaps help to explain the failure to 
detect any discrete cloud features from 45° S to 
the south pole. 

9. The failure of modelers of near-IR spectra to de¬ 
tect a methane cloud at the proper location (where 
the occultation found a sudden change is refractiv- 
ity) is due to two effects: one is that most mod¬ 
elers used mixing ratio profiles providing too lit¬ 
tle methane absorption and the other is that the 
methane cloud layer has a significantly lower opac¬ 
ity than the main cloud layer that peaks in the 
1.4-1.7 bar region, making it hard to resolve its 
separate contribution. 

10. The failure of the occultation measurements to de¬ 
tect the deeper and more significant middle tro¬ 
pospheric cloud layer that so prominently affects 
both near-IR and CCD spectra is not at all sur¬ 
prising if it is created by condensibles at very low 
mixing ratios. In that case the condensation would 
not produce a detectable change in refractivity. 

11 . Although it is not possible for the STIS spectra 
alone to clearly distinguish between vertically ex¬ 
tended models, such as that of KT2009, and com¬ 
pact layer models of the type we presented here, 
there is a strong distinction in physical processes 
involved. Our model is consistent with cloud for¬ 
mation by condensation, while the KT2009 model 
presents a picture in which a vertically extended 
haze of stratospheric origin provides all the particu¬ 
late opacity, much like the vertically extended haze 
on Titan. We prefer the condensation model be¬ 
cause it provides consistency with the occultation 
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observations and provides a mechanism for deplet¬ 
ing methane from the upper troposphere at high 
latitudes. 

Many of these results are dependent on models of 
methane absorption in spectral regions where absorp¬ 
tion is relatively weak and somewhat uncertain. Some 
of the absorption coefficients were adjusted by KT2009 
to obtain better overall fits to the observations, which 
might tend to favor the particular models they were us¬ 
ing in their analysis. As we obtain more accurate infor¬ 
mation about methane absorption that is independent of 
the spectral observations we are trying to interpret, these 
conclusions may need to be modified. 

The view of Uranus’ north polar regions at the be¬ 
ginning of 2011 is now adequate to address an impor¬ 
tant issue with new observations. The issue is whether 
the depletion of methane observed at southern high lati¬ 
tudes is also present at northern high latitudes, and if not 
currently depleted, whether it will become depleted as 
Uranus seasons proceed. There is some evidence suggest¬ 
ing high northern latitudes might not be depleted: dis¬ 
crete cloud features have been o bserved in 2007 Keck im - 
ages between 45° N and 74° N (|Sromovskv et al.ll2009D . 
while no discrete cloud features have been seen at cor¬ 
responding southern latitudes. The lack of southern dis¬ 
crete features is correlated with the inferred descent of 
dry gas at those latitudes, suggesting that the descending 
gas and resulting methane depletion may not be present 
at high northern latitudes. The question of whether 
methane is currently depleted at high northern latitudes 
on Uranus can be answered unequivocally by new STIS 
spectral measurements. 
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